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ABSTRACT 

We present a model for the formation of silicon carbide aggregates within the ex- 
panding and cooling supernova remnant. SiC-X grains measured in the laboratory at a 
high spatial resolution have been found to be aggregates of smaller crystals which are 
isotopically homogenous. The initial condensation of SiC in the ejecta occurs within 
a interior dense shell of material which is created by a reverse shock which rebounds 
from the core-envelope interface. A subsequent reverse shock accelerates the grains 
forward, but the gas drag from the ejecta on the rapidly moving particles limits their 
travel distance. By observing the effects of gas drag on the travel distance of grains, we 
propose that supernova grain aggregates form from material that condensed in a highly 
localized region, which satisfies the observational evidence of isotopic homogeneity in 
SiC-X grains. 

Subject headings: — supernova remnants — infrared:stars — astrochemistry 

1. Introduction 

An important source of interstellar dust is the rapidly expanding and cooling outflows of type 
II supernovae. Within the highly radioactive interior of the ejecta, Si, C and O atoms chemically 
react in an a nearly hydrogen-free environment, forming SiC, CO, and graphite. The most wel l- 



studied supernova condensate in the laboratory is the SiC type-X grains (IClavton &: Nittlerl 12 004) 



Recent studies of the structure and isotopic makeup of these grains (jStroud. Nittler. <fe Hoppdl2004l ; 



Hynes. Croat. Amari. Mertz. Bernatowiczll2006l ) show that SiC-X grains are typically comprised 



of an agregate of smaller isotopically homogenous subgrains a a few tenths of a micron to a micron 
in size. This microstructure indicates that dust coagulation processes likely occur within the ejecta. 



We follow iDeneault. Clayton. &; Hegerl (|2003l . hereafter DCH03) in proposing a model for the 



formation of SiC-X grains in the ejecta, in which SiC will condense in the dense ejecta interior, and 
be accelerated through the overlying ejecta via reverse shocks. SiC condenses within the dense, hot 
interior regions of the supernova ejecta. The reverse shock from the core-envelope interface sets up 
such a high-density region (between 2.7 - 3.6Mq). Figured] reproduces Fig. 2 from DCH03 to show 
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the large density enhancement which leads to SiC growth. Although the exact condensation kinetics 
of SiC within the hydrogen-free supern ova interior is currently unknown, the re is indication that 
strong shocks can drive SiC formation (jFrenklach. Carmer. &: Feigelson 1989 ). Within the dense 
region created by the reverse shock, SiC will only form when the Si/C ratio is greater than 10 
(DHC03), This condition in satisfied interior to 3.2M0 and we consider the region 2.7< m <3.2M 
to be the primary SiC-condensing region. 

The SiC grains' motion is coupled to that of the homologously expanding ejecta initially. A sec- 
ond reverse shock (between 10 7 - 10 9 s after core bounce) caused by the interaction between the fast 
moving ejecta and the surrounding circumstellar material causes these grains to be decoupled from 
the gas, and move forward through the overlying material. Recent models of grain survival in ejecta 
(Nath. Laskar, & Shull 2008; Noz awa. Kozasa. Habe. Dwek. Umeda. Tominaga. Maeda. &: Nomoto 



2007 



Bianchi &: Schneider 



20071 ) pose that upward of 20% of the dust mass will be destroyed by 



sputtering in the reverse shock. Smaller condensates are slowed by the drag force of the gas and 
non-thermal sputtering and are eventually destroyed by thermal sputtering while la rger condensates 
a few tenths of a micron in size survive roughly intact (IBianchi Schneiderll2007l ). 




Fig. 1. — Density of the ejecta in mass coordinate plotted for several different times. The large 
density enchancement between 2.7 - 3.6 Mq contains the SiC condensation region. (This figure 
replicates Fig. 2 from DCH03) 

Grains that condense farther out in mass coordinate will have a higher velocity than those 
which condensed closer to the center. Under this expectation, grains will typically not interact. 
However, DCH03 posed that larger grains formed in inner mass coordinates can overtake smaller 
grains that formed farther out, because the differential velocity of a grain due to gas drag is 
proportional to the surface area of the grain. By this method, larger grains would act to sweep 
up smaller grains in their path. The aggregation model of DCH03 ex pected that grains o f varying 
isotopic composition would aggregate together. Current evidence (e.g. iNittler et al.ll2007l ) suggests 
otherwise. Although the nanoscale, intragranular-resolution data on SiC-X grains is limited to a 
few grains, the available analyses suggest that grain aggregates are formed from sub-grains that 
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condensed from a single isotopic resevoir. 



2. Ballistic Motion Through the Ejecta 

2.1. Reverse Shocks 

The amount of pre-supernova mass loss determines the timescale for the reverse shock to enter 
the condensation region as well as the strength of the reverse shock. If the mass loss is too large, the 
shockwave will enter the condensation region while SiC growth is starting. The high temperatures 
in the shocked gas will stop condensation within the gas, and the small grains that may have already 
condensed will be destroyed by thermal and non-thermal sputtering. A reverse shock that arrives 
within the condensation region after SiC growth has slowed, t > 3 years, is more conducive to grain 
survival. Decoupled from the gas, larger grains are more apt to survive non-thermal sputtering, 



and are not melted by the shock (DCH03, iNozawa et al.1 120071 ). 



The inital grain velocity depends on where the grain condenses as well as the strength of the 
reverse shock. The initial grain velocity v 9 q is the difference between the homologous expansion of 
the gas and the velocity of the reverse shock. For the purposes of this paper, we take v g o to be a 
free parameter, and without loss of generality, we choose v g o = 500 km s _1 to be the "standard" 
initial velocity, following DHC03. 



2.2. Gas Drag 



Each grain is decelerated by a drag force which is caused by collisions between the grain and 
the overlying gas as well as non-thermal sputtering from collisions with gas particles. Non-thermal 
sputtering serves to erode the dust grain, slowing its forward motion. In this paper we neglect 
non-thermal sputtering, and thus provide upper limits to tr avel distances of grains based on gas 
drag alone. The grain velocity decreases due to gas drag by (jDraine Salpeterlll979l ): 
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where a g is the grain radius, and p g is the density of the SiC grain, and nj is the number density 
of species i. The function Gi(si) is summed over all gas species: 
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and Si is the atomic sp eed ratio Sj = (rriiV g /2kT) 1 / 2 . In this study, a s in DCH03, we use the 25M(7 







ejecta model s25 from iRauscher. Heger. Hoffman. &: Woosleyl ([2002J, , hereafter WHW). This is a 
one-dimensional ejecta model with solar metallicity, and omits instabilities (e.g. iHerant Woosley 
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1994 ) that can move inner-ejecta material to higher mass coordinates. The net effect of these 
instabilities on the ejecta is to move the condensation region for SiC grains to a different mass 
coordinate, but it does not significantly alter the motion of grains due to reverse shocks, or the 
condensation chemistry of SiC in these dense clouds. 

The effect of the gas drag on the motion of the grain through the overlying ejecta is pronounced. 
In Figure [2j we show the distance in mass coordinate traveled by the grain as a function of the grain 
radius and initial velocity through the ejecta. Here we have assumed the grain begins its motion 
at 3.OM0 in the s25 model. The total distance traveled is the distance that the grain has moved 
before the gas drag has slowed the grain enough to be comoving with the gas. Smaller grains are 
effectively and quickly slowed by the gas drag, slowing to a differential velocity of zero after only a 
few million seconds, effectlively trapping them within a localized region. Larger grains that formed 
interior to these grains will easily overtake them even though they are moving with a lower initial 
velocity. 




Grain Radius (/im) 



Fig. 2. — Distance in mass coordinate traveled by grains as a function of grain size for different 
initial velocities. Smaller grains, regardless of initial velocity remain near their local region, larger 
grains travel up to ten times further in mass coordinate. 



2.3. Grain-Grain Interactions 

In Figure O we see the path of a hypothetical 1 micron grain formed at 2.95Af . Three other 
grains that formed in the overlying layer at 2.99M© are also shown. For clarity, we have not shown 
grains that may have condensed between those two coordinates. The time axis has been scaled 
such that zero representes the time that the reverse shock reaches the condensation region of the 
grain. In the ejecta, the shock reaches the inner zone on the order of 10 s later. Although the 
interior grain has a slower initial velocity, its larger size minimizes the effect of gas drag. In the 
absence of sputtering or inelastic collisions, the larger grain can travel approximately O.IMq in 
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mass coordinate. The velocity of grain collisions decreases with increasing target size, because 
smaller grains are effectively smowed by gas drag. Low velocity (thus low kinetic energy) impacts 
have a much higher chance of grains stick ing and forming aggregat es compared with higher velocity 
impacts, which could shatter the grain. (jDominik &: TielensHl997l ). 




Fig. 3. — Larger grains formed interior to smaller grains can overtake them. The time axis has 
been scaled such that zero represents the time when the shock reached the grain. The grain that 
formed at 2. 95M & can interact with grains up to ~ 0. 50[im in radius that formed a few hundredths 
of a solar mass farther out. 

Although we have not included grain shattering in the calculation, it is worthwhile to consider 
the effect of grain shattering on the evolution of grains within the ejecta. The critical velocity for 
catastrophic destruction in a grain-grain collision is given by 

/i?|\ 9/16 

VcatCC [w P ) (3) 



Whe re Rt and Rp are the target and projectile radii respectively (adapted from lJones. Tielens. &; Hollenbach 



1996). In Figure [U we see that small projectiles, less than O.l^m need an extremely large relative 
velocity to shatter a target grain of any size. However, only larger grains can survive catastrophic 
spallation if the relative velocities between projectile and target are a few km s _1 . This implies a 
very stark picture for the survival of smaller grains. Larger SiC grains which condensed at inner 
mass coordinate will more likely shatter small grains unless the impact velocity is much less than 
a km s -1 . Aggregates of simil ar micron-size crystals have been observed in the laboratory (e.g. 



Nittler. Hoppe. fc Stroudll2007l ). These larger grains are more likely to survive impacts with other 
similar sized crystals at higher velocities. Fig. 0] shows that a l//m impactor will not catastrophi- 
cally shatter a 1/im target if the velocity is less than about 4 km s . 



Fig. 4. — Critical velocity for shattering a target grain as a function of radius of the impacting 
grain (Rt)- Five target radii are plotted from 0.1 to 1.0 fim. Larger SiC grains will shatter smaller 
grains at even low-speed impacts of a few tenths of a km s _1 . Only collisions between similar sized 
grains at low velocities will lead to aggregation in lieu of shattering. 



3. Condensate Size Distribution 



Absent any collisions with other grains, a micron sized SiC grain will travel not much more 
than 0.2Mq outward before gas drag slows its relative forward motion (compared to the ejecta) to 
zero. This distance represents a strict upper limit to the motion of the grain in the ejecta. If the 
condensation region for SiC grains is found between 2.7<m<3.2M0 (as discussed in JQ) that means 
that SiC grains can interact with material no further than m = 3.4M Q . In the inner region, near 
2.7M0, the Si/C ratio approaches 10 4 . In this highly Si-saturated region graphite will be readily 
converted to SiC (DCH03). Due to the much smaller carbon budget in this region however, it is 
very likely that the inner condens ation zone will produce both a smal ler number of condensates as 



well as smaller-sized condensates (jDeneault. Clayton. &: Meverl 12006). The specific kinematics of 



SiC formation as a function of the Si/C ratio are unknown, and we do not endeavor to solve this 
problem in the current paper. As a rule of thumb, we postulate that the size of SiC condensates 
increases roughly inversely proportional to Si/C ratio. 

As discussed in Q2.2\ the forward motion of small grains is strongly inhibited by gas drag. In the 
high-temperature post shock environment, thermal sputtering will efficiently destroy small grains. 
Since the interior condensation region produces only these smaller grains, aggregate formation is 
essential for SiC grain survival in this region. Because a smaller number of initial condensates are 
formed, however, aggregates from the interior region will be rare, and comprised primarily of smaller 
subgrains. Similarly, the Si/C ratio drops below 10 above 3.25Mq, which inhibits the growth of 
SiC. Grains that condensed near this outer edge of the region will not form SiC aggregates with 
overlying material, and will remain single crystals. The central region of the condensation zone, 
between 2.9 - 3.IM0, is the most likely place for aggregate formation. 



-7- 



4. Isotopic Composition of Grains 



DCH03 proposed that an aggregate of grains would have a noticable isotopic heterogeneity, in 
that each subgrain of the aggregate would have a markedly different isotopic co mposition. Labora- 



(2002); 


Nittler et al. 


( 


2007) 



mogeneity of supernova aggregates indicates that the process of coagulation occurs in a highly 
localized area of the ejecta. Due to gas drag, the travel distance of a grain is limited. It will only 
have the chance to interact and coagulate with grains that formed nearby. As discussed in §2.31 
large grains, which decelerated the least by gas drag will have only limited capability of interacting 
with overlying grains, smaller grains, much less so. 




Fig. 5. — Mass fraction of 2S,29,30 Si as a function of mass coordinate in the ejecta. We can see that 
the isotopic composition between 2.8 - S.IMq is fairly uniform. Below 2.8M & , the abundances of 
29 ' 30 Si drop by an order of magnitude. This demonstrates that grains that form near one another 
will have very similar isotopic composition. 

Does this interior region produce the 28 Si excess observed in SiC-X grains? As we can see in 
Figure [H the condensation and aggregation region between 2.9 - 3.1M does have a substantial 
28 Si excess. This excess varies by a factor of three over the region, but is always at least a factor 
of 10 higher than the abundances of 29,30 Sz. If we consider the supernova to be of half-solar 
metallicity (as per DCH03), the 28 Si abundance will be much higher in this region. The 14 iV/ 15 iV 
ratio, however, is near 1.6 x 10 -2 within the region. This is far lower than the measure d ratio in 
X grains, which sits between 20-200 ()Hoppe. Strebel. Eberhardt. Amari. Lewis boool ) . DCH03 
proposed that the implantation of 14 N would occur as the grain is propelled through the ejecta, but 
implantation is insufficient to overcome this problem, especially if the grain's travel is limited by 
the gas drag. Other isotopic ratios, such as 2e Al/ 27 Al, do not match the observed data. The 25M Q 
ejecta model that we have used in this paper does not consider the possibility of mixing of gases 
in the eiecta lYoshidal (|T.I . e.g.). Such ejecta mixing can help to explain the isotopic composition of 
the condensates. 
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5. Ejection into the ISM 

The rapidly decreasing grain velocities in our model effectively retard the total distance that 
grains can travel relative to the gas. This poses a severe problem, as the grains will remain effectively 
trapped within the expanding ejecta, and not ejected into the ISM. We propose that the grains will 
be ejected into the ISM by a third reverse shock between the outward moving ejecta and the ISM. 
The time for this reverse shock is much longer (10 9 - 10 11 s) than that of the previous shock due to 
the circumstellar material. At these late times, the ejecta is very diffuse, and the accelerated grains 
will not be slowed appreciably by the gas drag or by sputtering, and will emerge into the ISM. 

6. Conclusion 

In this paper we have presented a model for the aggregation of SiC grains in the supernova 
interior. In the one-dimensional model, SiC will condense within a limited region of the interior 
where the number density of Si atoms is greater than 10 times that of C. When the reverse shock 
passes through the condensation region, the grains are decelerated much less than the surrounding 
gas, and they ballistically move through the overlying ejecta. Gas drag from the ejecta slows these 
grains down with a rate inversely proportional to the radius of the grain. Larger grains can shatter 
smaller grains through impacts, but brains that have low relative velocities can aggregate together. 
The gas drag highly truncates the distance that a grain can travel, thus an aggregate will likely be 
comprised only of grains with similar isotopic composition. 

We have not set out in this work to show a detailed kinetic description of grain coagulation 
within the ejecta, but rather to provide a more refined physical model for that process beyond what 
was hypothesized in DCH03. The greatest uncertainty in this model is whether grains will actually 
coalesce in the low-speed ballistic collisions that we have presented, and what effect those collisions 
have on the evolution of the aggregate grain. We have also omitted the effects of sputtering and 
turbulent motion of the ejecta on the grains' travel. Future ongoing work in this area will provide 
us with a more complete picture of the processes of SiC dust aggregation within the ejecta. 
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